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Abstract 

The exact solutions of a one-dimensional mixture of spinor bosons and spinor fermions with 
^-function interactions are studied. Some new sets of Bethe ansatz equations are obtained by 
using the graded nest quantum inverse scattering method. Many interesting features appear in 
the system. For example, the wave function has the S U(2\2) supersymmetry. It is also found that 
the ground state of the system is partial polarized, where the fermions form a spin singlet state 
and the bosons are totally polarized. From the solution of Bethe ansatz equations, it is shown 
that all the momentum, spin and isospin rapidities at the ground state are real if the interactions 
between the particles are repulsive; while the fermions form two-particle bounded states and the 
bosons form one large bound state, which means the bosons condensed at the zero momentum 
point, if the interactions are attractive. The charge, spin and isospin excitations are discussed in 
detail. The thermodynamic Bethe ansatz equations are also derived and their solutions at some 
special cases are obtained analytically. 
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1. Introduction 

Recently, tremendous experimentalprogresses have been taken in the research of the one- 
dimensional (ID) trapped cold atoms |1 H TI H, EL S, S S S [HL E3] - By means of either 



magnetic or optical traps, the cold atom gas has been realized. With Feshbach resonance, the 
scattering length and thus the couplings among atoms can be manipulated. In addition, with 
laser beams, one can confine particles in the valleys of periodic potential of the optical lattice. 
These experimental tools provide a platform to study some controllable condensed matter sys- 
tems. Theoretical, many methods have been applies to study these problems. For example, the 
Gross-Pitaviskii (GP) theory is widely adopted in dealing with the systems of Bose-Einstein 
condensations with weak interactions. However, the GP theory is based on a mean field approx- 



imation and has many shortcomings. The boson-Hubbard model H 1 3L 1 1 411 are used to describe 
the low-energy physics of ultracold dilute gas of bosonic atoms in an optical lattice. This model 
can not be solved exactly and one has to use the numerical approaches or approximate methods 
such as the valence bond mean field lfl5ll . Exactly solvable models play an important role in 
the investigation of ID interacting many-particle systems. The exact solutions can supply some 
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believable results thus serve as a very good starting point to understand the new phenomena and 
new quantum states in trapped cold atomic systems. 

Using laser beams or Feshbach resonance, the local interactions can be strongly enhanced. 
The low-energy behavior of trapped cold atoms can be well described by a system with (5-function 
potentials. Fortunately, this kind of problems can be solved exactly. The exact solutions of scalar 
bosons with (5-function interactions are obtained by Lieb and Liniger 11611 . Li, Gu, Ying and 
Eckern study the two-component bosons, where the intrinsic degrees of freedom (isospin) of 
bosons satisfy the SU(2) symmetry They find that the ground state of the system is not 
an isospin singlet, but a polarized or ferromagnetic state. In fact, the Bethe ansatz equations 
of multi-component bosons have already been obtained by Zhou iflill . Naturally, one should 
consider the effects of spin exchanging interactions if the bosonic atoms have non-zero spin. An 
exactly solved model of bosons with spin-1 is proposed in the Ref. 11911 . where the exchange 
interactions are embodied in the model Hamiltonian and the system is spin-dependent. More 
physical properties of bosonic cold atomic systems with ^-function potentials can be found in 
theRefs. JaaElllllllllilll. 

In experiments, atoms with inner degrees of freedom (hyperfine spin) are prepared by catch- 
ing several hyperfine sublevels of atoms. The spins can be polarized by the external magnetic 
fields, thus it is easier to capture the inter state of the fermions by using the laser beams. Now, 
people can control the fermionic atoms steadily staying on some special inter states, and realize 
the fermionic atoms with multi-component hyperfine spin |9J]. The theoretical model of the spin- 
1/2 fermions with <S-function interactions is solved by Gaudin [26] and Yang [27]. Sutherland 
generalizes the results to the multi-component fermionic systems 12811 . 

Most recently, the study of ultracold Bose-Fermi mixtures become a remarkable research 
topic for many new matter phases may arise in these systems. Experimenters have succeeded in 
preparing the mixtures of 7 Li- 6 Li, 23 Na- 6 Li or 87 Rb- 40 K in the optical lattices ll29ll30l l3ll 1321 133, 
34,135, 3a 37]. For example, a stable bosonic 87 Rb and fermionic 40 K mixture in three dimen- 
sional optical lattices has been realized Bo, 13711 . Theoretically, Lai and Yang obtain the exact 
solutions of ID mixture of spin- 1/2 fermions and scalar bosons with the 5-function interactions 
JHI]. They calculate the ground state energy and gapless fermionic excitations in the thermody- 
namic limit 13911 For more studies on the boson-fermion mixtures in the optical lattices, please 
see the Refs. J4flEl 53 • 

It is natural that one should consider the mixtures of multi-component fermions and multi- 
component bosons. In this case, both the fermions and the bosons have the intrinsic degrees of 
freedom. It is well-known that the ground state of bosonic systems with intrinsic degrees of free- 
dom can be surprisingly different from that of the scalar bosons. Furthermore, the ground state of 
the spin- 1/2 fermionic system with ^-function interactions is spin singlet, while the ground state 
of bosonic system is isospin polarized or ferromagnetic state. One may wonder what will happen 
if we mix the bosons and fermions both with some intrinsic degrees of freedom? What is the 
new quantum state in the mixtures? These issues are quite interesting and important nowadays 
due to the rapid progress in the field of cold atomic physics. 

In this paper, we study a mixture of two-component fermions and two-component bosons 
with ^-function interactions. Because the wave function of the system is supersymmetric and 
satisfies the S U(2 12) Lie sup eralgebra, we use the super or graded nest quantum inverse scattering 
method |44, 45, 4a, 47, 48, 49[ 50*] to derive the exact solutions of the system at ID. We obtain 
the Bethe ansatz equations with different gradings. We find that the ground state of the system 
is partial polarized. If the interactions are repulsive (c > 0), all the momentum, spin and isospin 
rapidities in the solutions of Bethe ansatz equations are real. If the interactions are attractive 
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(c < 0), the Bethe ansatz equations may have the complex solutions, where the fermions form 
the two-particle bounded states and the bosons are condensed at the zero momentum point at the 
ground state. We then study the low-lying excitations such as charge, spin and isospin excitations 
in detail. We also obtain the thermodynamic Bethe ansatz equations at finite temperatures and 
find their analytic solutions at some special limiting cases. 

The integrable S t/(2|2)-supersymmetry is a very interesting issue. The Bethe ansatz for 
the corresponding quantum spin chain is obtained by Kulish 14411 . and the continue limit and 
low energy behaviors are studied by Saleur II50I1 . In this paper, we study the corresponding 
continue quantum gas model, which also has many applications in the systems of cold atoms 
with hyperfine structure. The low-lying excitation spectrum can be measured by the scattering of 
phonons in experiments. Meanwhile, the cold atoms with different internal states can be prepared 
in experiments. T he p honons scattering experiments on the bosonic cold atoms have been done 



several years ago 15 1L 15211 ■ Thus the motivation of this paper is to give a prediction on the ground 
state and low-lying excitation properties of the bose-fermi mixture in the cold atom systems, for 
the excitations in spin and charge sectors can show different behaviors. 

The paper is organized as follows. In section 2, we introduce the supersymmetry of the 
system and the generators of the corresponding Lie superalgebra. In section 3, we derive the exact 
solutions of the system by using the generalized quantum inverse scattering method. We give the 
Bethe ansatz equations with the BBFF grading, which are the foundations of our discussions. The 
ground state properties are discussed in section 4 and the low-lying excitations are discussed in 
section 5. The thermodynamic Bethe ansatz equations are calculated in section 6 and some useful 
limit cases are discussed in section 7. Section 8 contains some conclusions and discussions. 



2. Supersymmetry of the system 

We study a ID cold atomic system mixed by bosons of species 1, Nb 2 bosons of species 
2, Nj { fermions of species 1 and N/ 2 fermions of species 2. The Hamiltonian reads 

XL r-L 
^ J d x cf>l(x)d x <pa(x)dx+ I ^g a/3 (pl(x)^ fj (x)ct) l3 (x)^a(x)dx, (1) 
a J° af> 

where a,/3 = b\,bi,f\,fi and <p a are the bosonic or fermionic field operators. The bosonic field 
operators satisfy the commutation relations, [<pl(x), fyiy)] = 5 a p5 xy , while the fermionic field 
operators satisfy the anticommutation relations, {(pl(x), <pp{y)} = 5 a p5 xy . According to the Pauli 
exclusion principle, the a and j3 in Hamiltonian (Q~|) can not be the same species of fermions. 
In this paper, we use the periodic boundary conditions. The wave function of the system (JTJ 
is supersymmetric, *P(x;, x/) = Pji¥(xi, xj), where Pfl means exchanging both the coordinates 
and the spins (isospins) of two particles j and /. Pji = 1 for bosons or bosons and fermions 
thus the wave function is symmetric, while Pji = -1 for fermions thus the wave function is 
antisymmetric. 

The supersymmetry of the system (Q]) can be described by the S U{2\2) Lie superalgebra. The 
superalgebra S U(2\2) has fifteen generators and eight of them are fermionic 



Moreover, the pure two-component fermionic subsystem has the S U(2) invariance and the pure 
two-component bosons also has the S U (2) invariance. In order to explain the symmetry of the 
system more clearly, we introduce the particle creation (annihilation) operators as a^-(x) (aa-(xj), 
where we assume two species of fermions carrying the different spins cr =[, |. There are four 



3 



kinds of states at a given position x, 

|0>*, I T)x = fl}(*)|0>„ \ i) x = a\(x)\0)x, \ / [i) x = a\{x)a\(x)\Q) x . (2) 

The state |0) x is vacuum and the state | |X) represents that an atom-pair is localized on a single 
energy level. Now, we introduce the generators of the superalgebra S U(2\2). The spin operators 
are defined as 



' + = J" a^(x)ai(x)dx, S = J" aj 



(x)a-i(x)dx, 



1 



■L 



S z = — I [al(x)at(x) — a. (x)ai(x)]dx. (3) 
2 Jo 1 1 

The spin operators S + , S ~ and S : form the S U (2) Lie algebra, where the commutation relations 
between the generators are [5 ,_ ,5' + ] = 2S Z ,[S + ,S Z ] - S + ,[S~,S Z ] = —S~. The above spin 
operators are grassmann even (bosonic). We introduce the pairing operators, which are also 
bosonic generators, 

if — I a] (x)al(x)dx, i] — I ai(x)aj(x)dx, 
Jo 1 1 Jo 

, 1 r L t t 

rf = - J [1 - fl T (x)a T (x) - a i (x)a i (x)]t/x. (4) 

The pairing operators rf,rf, rf also form a S U(2) Lie algebra with the commutation relations 
[t]~, r/ + ] = 2rf, [t] + , rf] = r] + , [rj~, rf] = —rf. The eight fermionic generators are 

Q<r = f [1 - at (x)ag-(x)]ao-(x), Ql= I [1 - at (x)ag-(x)]al.(x), 
Jo Jo 

Q,r = J at (x)aa-(x)aa-(x), Ql = I at(x)ag-(x)al(x), (5) 

Jo Jo 

where cr means the spin with opposite direction of <x. The fermionic operators are grassmann 
odd. These operators together with the unit operator J Q Idx = L generate the Lie superalgebra 
SU(2\2). 



3. Bethe ansatz solutions of the system 

In the following, we consider the case that all the coupling parameters are equal, g a p = c. 
The system ([TJ has several integrable lines. If Nt-, = N/ t = Nf, = 0, the system degenerates 
to the scalar bosons with (^-function interactions which is solved by Lieb and Liniger 01611 . If 
Nj\ = Nf = 0, the system degenerates to the two-component S U(2) bosons and is studied by Li, 
Gu, Ying and Eckern 111 711 . If = = 0, the system degenerates to the spin- 1/2 fermions 
and is solved by Yang II27I1 . If = 0, the system degenerates to the mixture of scalar bosons 
and spin-1/2 fermions which is solved by Lai and Yang [ 38] . In this paper, we consider the case 
that all the particles numbers , Nb 2 , Nfi and N/ 2 are not zero. We fist derive the two-body 
scattering matrix by using the coordinate Bethe ansatz method and prove the integrability of the 
system. Then we determine the Bethe ansatz equations and the energy spectrum by using the 
nest quantum inverse scattering methods. 
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3.1. Coordinate Bethe ansatz 

In the framework of coordinate Bethe ansatzjhe wave function of the system described by a 



set of quasi-momenta {kj} can be written as B271I38I1 



••• ,x N s N ) = J^GixQ, < ■ ■< x QN )A Sl ... SN (Q,P)e i ^ k n x e l , (6) 
Q.P 

where Q = (gi, ••■ , Qn) and P = {P\, ■ ■ ■ ,Pn) are the permutations of the integers 1, ■ • • ,N, 
N is the total number of particles, N - + Nb n + N/ t + N/ 2 , 6(x Ql < ■■■ < x Qhl ) = 6(x Qn - 
x Qn-i ) ' ' ' 8( X Q^ ~ x Qi ) an d 6( x ~y) is the ste P function. The wave function is supersymmetric under 
permutating both the coordinates and the spins (or isospins) of two particles. The wave function 
is continuous but its derivative jumps when two atoms touch. With the standard coordinate Bethe 
ansatz procedure, we obtain the two-body scattering matrix as 

kj -k,- icP s . 
kj + (?) 

where Pj { is the spin super permutation operator with the definition [P^]^ = i~^) e ^ ep S av 6 ll p, 
the a and \x are the row indices, and B and v are the column indices. Here e a is the grassmann 
number, e a — for bosons and_& = 1 for fermions. The scattering matrix satisfies the super or 
graded Yang-Baxter equation lHl^|46ll4^ l48ll49l l50ll 



S n(ki - k 2 )S n(ki - k 3 )S 23(^2 - h) = S 2 3(ki - k 3 )S l3 (k\ - k 3 )S 12(^1 - k 2 ). (8) 

which ensures the integrability of the model (HJ. The Yang-Baxter equation (O can also be 
written out explicitly as 

s - k 2 t\%s 13 (t, - k 3 y b \%s 23 (h - k 3 ) c b 2(-)^ + ^- 

= S 13 {k 2 - k 3 ) b £]S 13 (*i - ^3)^^12(^1 - £ 2 )^ 2 2 (-) (£ "' +e <" ^ . (9) 

With the periodic boundary conditions of the wave function, we obtain the following eigen- 
value equations 

S jN(kj - kw)S jN-\{kj - fev-i) ' ' ' S jj+i(kj - kj+\) 

xSjj-^kj - kj-i) ■ ■ -Sflikj - kOe ik J% = ft, (10) 

where £0 is the amplitude of initial state wave function. 

3.2. Algebraic Bethe ansatz 

Now, we derive the exact solutions of the system by using the graded nested quantum inverse 
scattering method 1 44 , 45 , 46[ 47[ 5H H2I ■ We consider the exact solutions of the system with 



BBFF grading, that is the grassmann parities for the four bases are e\ — e 2 — and e 3 — 64 = 1 . 
Please note choosing different bases is equivalent to choosing different highest weight represents 
when deriving the Bethe ansatz equations. The Bethe ansatz equations with different gradings 
can change into each others by using some transformations |48]. 
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The matrix form of the scattering matrix Sqj(A) in the space is 



a(A) - b(A)e) 



-b(X)e 
-b(A)e 
-b(A)e 



-b{A)ef 
a(A) - b(A)ef 
-b(A)ef 
-b(A)ef 



-b(A)e 
-b{A)e 



a(A) + b(A)ef 
b(A)ef 



-b{A)e 
-b(A)e* 
b(A)ef 



.41 



a(A) + b(A)e*j 



44 



(ID 



where the matrix e"/ 1 acts on the y'-th space with its elements defined as (e"^) MV = 5 afI 5p v , a(A) = 
A/(A + ic) and b(A) = ic/(A + ic). The quantity (fTTT) is also called as the Lax operator acting on the 
j-th space. Introduced the braid scattering matrix R\ 2 (A) = P s n S 12(A), which satisfies the braid 
Yang-Baxter equation, 



R n (A - u)R 23 (A)R l2 (u) = R 23 (u)R l2 (A)R 23 (A - u). 



(12) 



We follow the graded nested algebraic Bethe ansatz method to solve the eigenvalue equation 
( [Tol l. The monodromy matrix is defined as 



T N (A) = S ( )j(A - kj)S N(A 

An (A) A 12 (A) 

A 2l (A) A 22 (A) 

A 3 i(X) A 32 (A) 

CM) C 2 (A) 



- kff) ■ ■ -Soj+\(A - kj +1 )S oj-\(A ■ 

A 13 (A) Bi(A) ^ 

A 23 (A) B 2 (A) 

A 33 (A) B 3 (A) 

C 3 (A) D(A) 



kj-i)---Soi(A-ki) 



(13) 



where means the auxiliary space and I = 1,2, ■ • -N mean the quantum spaces. From Eq. (1121) . 
we can prove that the monodromy matrix ( fT3l satisfies the Yang-Baxter relation 



R l2 (A - u)[T N (A) ®s T N (u)] = [T N (u) ® s T N (A)]R l2 (A - u), 



(14) 



where ® s means the super or graded tensor-product as [A <g>. s B~\ h a d c = (-\)^" +eb)<Lc A a i,B cc i, a and c 
are the row indices, and b and d are the column indices. Using indices, the Yang-Baxter relation 
<TT~4T > can also be written as 



R n (A-u) b a ^T N (A)l[T N (u)f 2 (-l) 

= T N (u) b a \T N (X)%R l2 {A- U )l\l(-iy^ K -- 



(e b , +ec, )e bl 



(15) 



\bib 2 



where all the repeated indices should be summed. The elements of scattering matrix S ij(A) a \ ai 
are not zero only with the conditions (1) a\ — a 2 — b\ — b 2 or (2) ct\ = b\,a 2 = b 2 or (3) 
a\ = b 2 ,a 2 = b\. These properties will be used in deriving the commutation relations. The 
transfer matrix t(X) of the system is defined as the supertrace of the monodromy matrix ( fT3l ) in 
the auxiliary space, 



t(A) = strT N {A) = Yj.-ir-TNiAt = A n (A) + A 22 (A) - A 33 (A) - D n (A). 



(16) 



From the Yang-Baxter relation (TPfl i. we can prove that the transfer matrices with different spectral 
parameters commute with each other [t(u), f(v)] = 0. Thus the system has infinite conserved 
quantities and is integrable. The eigenvalue problem (fTOl is therefore reduced to 



- str Q T N (kj)e ik >% = ft. 
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(17) 



We choose the local vacuum state as |0) ; - = (0,0,0, 1)' where t means the transpose. The 
global vacuum state is constructed as |0) = ® N j =i \Q)j. Acting the monodromy matrix ( TT3T > on this 
vacuum state, we have 



T N (A)\0) = 








^ 

Uti a(A -h) 

\\t\ a{ - A - k i) 

C 2 (A) C 3 (A) 1 



10). 



(18) 



We see that the elements An (/I), A2i{A), A^{A) and D(A) acting on this vacuum state give the 
eigenvalues. The elements B a (A) acting on the vacuum state are zero. The elements C a (A) acting 
on the vacuum state give nonzero values and can be regarded as the creation operators. We 
assume the eigenstates of the system (Q3 are obtained by applying the creation operators C a on 
the vacuum state as 



Hi, • ■ • , A Nl \F) = C ai (Ai) ■ ■ ■ C am (A Nl )|0>F fl "> "*> , 



(19) 



where F a "' is a function of the spectral parameters Aj and Ni is the number of creation opera- 
tors. When the transfer matrix acting on the Bethe states (I191 l. we need the commutation relations 
between An, A22, A33, D and C a . From the Yang-Baxter relation ( fTSt and using the properties of 
the R matrix, we find following commutation relations 

1 . b(A-u) 



D{u)C c (A) 



-C c (A)D(u) ■ 



a(A - u) 
-(-])' 



a(u - A) 

l f „4-l)fe+l)^("-^) 



( -C c (u)D(A), 
a(A — u) 

C e {A)A ad (u) 
C h (u)A ac (A), 



a(u - A) 

C ai (u)C a2 (A) = R%b( u ~ ^C t! WG,(4 
where all the indices take values 1, 2 and 3. The first nesting R matrices are defined as 



R ( p FB (u) = b(u) + a(u)P [ 



m 

FFB' 



-b(u) + a(u)P [ 



,(i) 
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where P FFB and P BBF are the 9x9 super permutation matrices for the grading e\ 



(20) 

(21) 
(22) 

(23) 




ex = 1, £3 

and e\ — €2 — 0, £3 = 1 , respectively. 

Acting the transfer matrix dT7b on the assumed eigenstate (T% , applying repeatedly the com 
mutation relations d20b - d22"l > and using the result ( fT8l , we have 

In, . n 

t(u)\Ai , ■ ■■ A Nl \F) = I Ff — Ff a(u - kt)t m (u) 

15=1 a(u ~ Aj) '=1 



-n 



f a(Aj — u) 



\Au--- ,A N] \F) + u.t., 



(24) 



where t (l \u) is the first nesting transfer matrix and u.t. means the unwanted terms. If the un- 
wanted terms cancel with each other, which gives following Bethe ansatz equations 



n 



n 



1 



a(A a - Aj) 

a(Aj - A a ) 1 = 1 a(A a - k t ) 



AD 



M-b 



(A V 1 F aN i 



1,2, 



(25) 



then the assumed states (TT9b are the eigenstates of the transfer matrix t{u) and the corresponding 
eigenvalus are given by the first term in Eq. (l24l >. 

Now, seeking the eigenvalues of t(u) becomes seeking the eigenvalues of t (l \u). The elements 
of first nested transfer matrix ? (1) (m) can be written out explicitly 



b v -b N } ^ Nl a N ( i) 



(26) 



Here all the indices c, are summed over and = P B l BF R B l BF {u). ^ n or der to interpret f l) {u) 

as the supertrace of the monodromy matrix, we define a new graded tensor product [F®G] h a d c = 
F*G^(-l) <£o+Ei,)(e ' +1) . This new graded tensor-product switches even and odd grassmann parities. 
Meanwhile, the first nesting monodromy matrix is defined as 



T^(u)=S^(u-A, 



' 01 



(u-Ai) 



' Ajjjdi) Ag(«) B ( }\u) ' 
= A^(u) Af(u) Bf{u) , 
, d/V) C™(w) , 

which satisfies the graded Yang-Baxter relation, 

Hu - v) [TV(u)®T$(v)] = [T«(v)®r™( M )] r(« - v), 

where the r-matrix is /*? = -b(u)6 a b5 c d + a(w)(W^c(-l) £ " +et+e " e ' ■ Then the transfer matrix 
is the supertrace of the first nesting monodromy matrix d27T i 



(27) 



(28) 



b,--h N 



strT ( V(u) = A ( V(u) + A<V(u) - D (1 \u). 



(29) 



We choose |0) 



(i) 



(0,0, 1)' as the local reference state for the first nesting. The global 



reference state is |0) (1) = ® ! Jl l \0) ( j > - Acting the first nesting monodromy matrix (l27T i on this 
reference state, we have 











o nfi < u - A i) 

Cf\u) C ( 2 l) (u) 1 J 



!0> 



(i) 



Assume the eigenstates of the first nesting transfer matrix are 



(1) 



b N ^N 2 ' 



(30) 



(31) 



where N2 is the number of creation operators. From the Yang-Baxter relation ([28), we obtain the 
following commutation relations 



D m (u)C^(A) 



a(A - u) 



a(A — u) 



(32) 



A^CfiA) = R °t~ A f e C?\A)A% (») + ^^Ci 1 >( m) A£(A (33) 



2< 1 )(M - /I) 



a(w — A) 



C^(u)C ( b \\A) = J?® (ii - A) b £c%(A)CW(u). 



(34) 



(2) 

Here, all the indices take values 1 and 2. The second nesting R matrices are R FF (u) — b(u) + 
a{u)Pp F and R^ B (u) = -b(u) + a(u)Pg B , where P FF and P BB are the 4x4 super permutation 
matrices for the grading e\ = ex = 1 and e\ = €2 = 0, respectively. [P^p]^ - -5 a d5bc, IP^ac = 
5 a d&bc- Acting the first nesting transfer matrix (1291 on the assumed states PIT) , we have 

"' |G)= n^^n a(M_ " /)f(2)(M) 

[ ;=i a(.u - Aj ) /=1 

"A -TOT- }l^ 1) .-.4 ) |G> + ^ (35) 

where f (2) (t<) is the second nesting transfer matrix 

? (2) («) = strSfl(u - A^Sfl^u - A^) ■ ■ -Sfl(u - a\\ (36) 
If the unwanted terms cancel with each other, which give following Bethe ansatz equations 

FT —L Li FT __J G W>» = PKA^t^G ^, a=l,2,---,N 2 , (37) 

the assumed states OTb are the eigenstates of the first nesting transfer matrix ( f29T >. Now, the 
remanent problem is finding the eigenvalues of second nesting transfer matrix ft\u). 

The second nesting scattering matrix is S {1 \u) = a{u) - b(u)P BB . The second nesting mon- 
odromy matrix is 

t { V<„\ - cG9 iOKcP) /-„ ^ ?( 2 V„ j^h - / A<2) ( M ) # (2) ( M ) \ n o~. 

i A , 2 (M)-i 0A , 2 (M-/i A ,Ji (W2 _ 1 (M-/i A , 2 _ 1 )---i 01 (K-/t 1 )-| C (2) (m) D (2)( m ) I. (3«) 

which satisfies the graded Yang-Baxter relations 

- v)[Tg(ii)®r®(v)] = [T$(v)®T™ (u)]R^(u - v), (39) 

where R (2 \u) = P^ BB S (2) (u) is the second nesting braid matrix. Here the graded tensor-product 
degenerates to the ordinal tensor-product because that only the bosonic gradings are left. The 
second nesting transfer matrix is the supertrace of the corresponding monodromy matrix 

f (2) («) = strT®(u) = A (2) (u) + D (2 \u). (40) 

We choose |0)^ 2) = (0, 1)' as the local vacuum state for the second nesting. The global vacuum 

state is |0) (2) = ®^j|0)* 2) . Acting the second nesting momodromy matrix d38l l on this vacuum 
state, we have 

T N2 W\0) -\ C(2)(M) U^a^u-A^-b^u-Af^r ■ (41) 

Assume the eigenstates of the second nesting transfer matrix f 2 \u) are 

\Af\ ■■■,41) = C^\A?) ■ ■ ■ C^)|0/ 2 >, (42) 



where Nj is the number of creation operators. From the Yang-Baxter relation d39b . we obtain 
following commutation relations 

D®(«)C®tf) = a{A - U) - b{A - U) cV\A)DV\u) + b ^-^C^\u)D^\A\ (43) 
a(A — u) a(A — u) 

A^(u)C^(A) = «(«-^)-M«-^) c(2)wa(2)(m) + b^^^ m 
a(u - A) a(u - A) 

C (2 \u)C (2) (A) = C (2) (A)C (2 \u). (45) 
The second nesting transfer matrix f 2 \u) acting the assumed states ( l42b gives 

[ N * a(u - A (2) ) - b(u - A (2) ) N 2 

P\u)C*\A?) ■ ■ ■ C«(jg)P>« = [1 ; () J n a(u - A?) 

[j=i a(u-A\>) )J 

*a(Af-u)-b(Af-u)* ) ... 

+ n t& i — n «<« - - b(u - ^} ^ ■ ■ • • • % i> + uj ~ (46) 

7=1 a \ A j /=1 j 

If the unwanted terms cancel, the assumed states (l42l are the eigenstates of the second nesting 
transfer matrix ? (2) (m), which gives following Bethe ansatz relations 

* a(Af> - Af ) - b(A {2) - Af ) a(Af> - J< 2) ) _ * fl(<l £> - jW) - fc (/t g> - ^) 
JJ ff a(Af - Af) - b{Af - Af) a(Af - Af) = \j a(Af - Af) ' m 

where a — 1 , 2, ■ • • , JV3 . Now, the eigenvalues of the transfer matrix f (m) can be calculated directly 
by synthetically considering Eqs. d24l i. d35l l and d46b . 

The forth set of Bethe ansatz equations are obtained from the eigen-equation ( TTOb as 

N, 

e- [k > h = \\ a(A a - kj), j = 1, ■ • • , 7V, (48) 

Put Aj — > /lj - /c/2, Aj — > /l^ - z'c and — > /i® - ic/2, then the Bethe ansatz equations 
(l25i (|37] >. (|47] | and d48]l can be written out explicitly as 



= 1 ""J 

Aa-*l-£c ^ A, a -A p -ic^Af-A a -{c 

1 jO) _ 3 777' 

y=l 7 a 2 



1 = 1 7^7-777 = 11 A a -A p + ic 11 



a =!,•■• ,7Vi, (50) 



.a) _i 777 = U 1W-W-T ' P = h-,N 2 , (51) 



7=1 A p A y + 2 L p=l A P *)S 2° 



^ 4 2) - 4° + ic J-i, A (2) -A (2) -ic 

r u A) , = n ; 2) 5r- ? r=i.-.^. (52) 



f=i A y A f 2 C i?=i,*r ^ y 
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where N = N bl + N bl + N fl + N fl , N\ = N b , + N h + N fi , N 2 = N hi + N b2 , N 3 = N b , . Taking the 
logarithm of Eqs. d49b - (p2l l. we arrive at 



kjL = 2nlj - J] E 1J2 (kj - A b ), j=l,---,N, (53) 

b=l 

N N, N 2 

2nJ a = Yj S i/2W« S i(^ - ^) + 2 El/2{Aa ~ fl = *» ' ' ' ,JVl ' (54) 

/=1 b=l c=\ 

N, N 3 

2^; ) = ^H I/2 (4 1 i ) -^)-^]H I/2 (4 I i ) -4 2) ), fll = l,...,^ 2 , (55) 

Z>=1 c=l 

2nJ% = - £ a 1/2 (4« - 4") + I] S l( 4 2 > - 4 2) ), «2 = 1, • • • ,N>, (56) 

6=1 c=l 

where S m (jc) = 2tarT 1 [x/(mc)], 7y, and are integer or half-odd quantum numbers. 

Here we have used the formula lnlYfc + ic)/(k - ic)] = i[n -2taiT l (k/c)]. These equations are the 
special case of that in the Ref. Il 811 . If we set N bl = N bl = 0, the system (1) degenerates to the 
spin- 1/2 fermions model with (5-function potentials and our Bethe ansatz equations d49b - d52l 



are the same as the ones obtained by Yang B27I1 . If we set N bl = 0, the system (1) degenerates to 



the bose-fermi mixture considered by Lai and Yang IB8I1 . and our Bethe ansatz equations are the 
same as their results. 

The eigenvalues of the Hamiltonian (1) are E = kj> where possible values of the mo- 
mentum kj are determined by the Bethe ansatz equations (1531 - ( |56T ). 

3.3. Corresponding lattice model 

From the transfer matrix t(u), we can construct the corresponding lattice model by using the 
standard integrable theory in the statistic physics. The Hamiltonian can be obtained by taking 
the derivative of the logarithmic form of the transfer matrix at the zero spectral parameter point 



05 311 . After some calculations, we find that the d\nt(u)/(du)\ u= o gives the Hamiltonian studied by 



Essler, Korepin and Schoutens 11461 l47l l48l 14911 up to a constant. The system (HJ is the continue 
case of Essler- Korepin-Schoutens (EKS) model. The two-body scattering matrix (fTTT > of present 
model has the same structure as that of EKS model. In such a sense, the spin dynamics of present 
model keeps some similarities to that of EKS model. 

We note that the integrable S U («|m)-supersymmetric quantum spin chain has been studied by 
Kulish [44]. The corresponding Hamiltonian can be obtained directly from the transfer matrix. 
If the charge rapidities tend to zero, our results degenerate into that obtained by Kulish. The 
thermodynamics and low-energy limit of the S f/(n|m)-invariant spin chain are studied by Saleur 
lf30h . These results are also valid for the present model in the spin sector. 



4. Ground state of the system 

We first consider the case of c > 0, which means the interactions among the particles are 
repulsive. We use the BBFF grading as a demonstration. Above we have confined the particles 
in a finite ID box with the length L. Some useful properties can be obtained by the analysis of 
the poles or zeros of the Bethe ansatz equations in the thermodynamic limit, that is the system 
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size L, particles numbers N,N\,N 2 and N3 tend to infinity, but the ratios N/L,Ni/L,N2/L and 
N3/L keep finite. For example, if some kj are in the upper complex plane, then the left hand 
side of Eq. d49l tends to zero when the system size tends to infinity. Thus the right hand side 
of Eq. (|49b should go to zero too. From further analysis of the Bethe ansatz Eqs. d49t - ( f52l l. 
we find that the momentum kj, rapidities A and A (V) are real at the ground state. The particle 
number = thus the corresponding rapidity A^ is zero too. That is to say only one species 
of bosons are left or the bosons are totally polarized at the ground state. The ground state for 
the bosons is a ferromagnetic state. Meanwhile, the particle numbers of two species of fermions 
are equal, N/ t = N/ 2 and they form the spin singlet states. The ground state for the fermions is a 
antiferromagnetic state. Therefore, the ground state of the system (1) is partial polarized. 

The momentum kj, rapidities A and A m at the ground state satisfy the following coupled 
equations 

N 



k j L = 2nIj-Y J Zv2(kj-A b ), (57) 

6=1 

N N, N 2 

2nJ a = ^ Si/aft, S i^« - A b) + J] S i/2(^ - (58) 

1=1 6=1 c=l 

Ni 

2nJ^ = Yj R y^ ~ h). (59) 



The quantum numbers Ij take integer (half-odd integer) values if Nb, + Nf, is even (odd), J a take 
integer (half-odd integer) values if Nb, + N/ 2 + 1 is even (odd) and take integer (half-odd 
integer) values if Nb, + Nf, is even (odd). Then quantum number configuration of the ground 
state are {/,■} = {-(N - l)/2, -(N - 3)/2, ...,(N- l)/2}, {/„} = {-(AT, - l)/2, -(N\ - 3)/2, . . . , 
(Ni - l)/2} and [J™] = {-(N 2 - l)/2, -(N 2 - 3)/2, ...,(N 2 - l)/2}, which are symmetrically 
centered around the origin if Nf, = N/ 2 is odd and Nb = Nb, is even. 

In the thermodynamic limit, the summations become integrations. The quantum number 
Ij, J a and become continue functions of the spectral parameters k, A and A m , respectively. 
Denote the densities of momentum k, rapidities A and by p(k), p\{A) and p^{A m ), respec- 
tively. Then we have p(k) = dIj/(Ldk), pi(A) = dJJ(LdA) andp (1) (/l (1) ) = dJ ( a \ } /(LdA m ). Taking 
the derivative of Eqs. (57\ - j59l , we obtain the densities of states at the ground state as 

m = — I f B 2c P^ dA 
Pi ) 2n + n J_ B c 2 +4(k-A) 2 ' 

1 r Q 2cp(k)dk 1 r B cpi(A')dA' 1 C D 2cp m (A m )dA {V) 

- - J q c 2 ~ 4(/l _ fe) 2 _ ~ J_ b c 2 + p _ ^)2 + ~ J_ D c 2 +4( / l_ / l(D)2' 



F 7rJ_ B c 2 +4(, 



i(A)dA 



■ 4(A^ - A) 2 ' 
The integral limits Q, B and D are determined by 



(60) 



rQ rB rD 

N/L= p(k)dk, Ni/L= pi(A)dA, N 2 /L= p m (A m )dA m . 

J-Q J-B J-D 



The quantity Q is the Fermi surface of the system. The densities of energy and momentum at the 
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ground state are 

rQ rQ 
E/L= k 2 p(k)dk, P/L= kp(k)dk. 

J-Q J-Q 

The magnetization of the fermions Sy and the magnetization of bosons T z are 
S z f = N fi - N f2 = 0, = N bl - N bl = N bl . 

In the case of c < 0, besides real solutions, Eqs. d49l ) - d52b also have complex solutions 
which are usually called as string solutions. After some algebraic calculations, we find that at 
the ground state, the momentum kj have the following 2-string and N bl -string solutions, 

kj = Aj+ l ^(3-2j)+o(e- 6L ), 7=1,2, (61) 

k, = ^(N bl +l-2l) + o(e- s ' L ), 1= 1,2 ■■■,N bl , (62) 

where A ; is a real parameter, 6 and 6' are some positive constants. Eqs. doTt and d62l means that 
the fermions form the spin singlet states and the bosons condensed at the zero momentum point. 



5. Low-lying excitation 



In this section, we consider the low-lying excitations in the system. We will follow the 
methods proposed by Takahashi 154 D5i 15611 . Essler and Korepin [49] very closely. The low- 
lying excitations are very rich due to the complicated solutions of the Bethe ansatz equations. 
Let us consider them one by one. 



5.1. Charge-hole excitation 

The simplest excitation is obtained by removing a quantum number Ij from the sequence {Ij} 
and putting it outside the sequence, i.e., 

{Ij} = ' m ~ hm +!,■■■ ,^y-!-,/„j, 

where /„ = (N -l)/2 + n, and keep the other two quantum number sequence \ J a , } unchanged. 
We call this excitation the charge-hole excitation since a "hole" is created under the Fermi surface 
and a particle outside the surface. The charge-hole excitation spectra are shown in Fig. [T] The 
dispersion relations of charge and hole are shown in Fig. [2] 

In the thermodynamic limit, the densities of states at this excited state read 



p(k) 



= J_ + 1 f* 2c Pl (A )d A 1 

2n n J_ B c 1 + 4(k - A) 2 L V ; 



1 r u 2c P (k)dk 1 r B cp x (A')dA' 

~ n J-Q c 2 +4(A-k) 2 ~ n J_ B c 2 + (A - A') 2 



1 r u 2c 

X J-D C 2 

1 r B 
W = - 

n J-B 



D 2cpV(AV)dA m 1 2c 



+ 4(A - /I* 1 )) 2 nL c 2 + 4(A - kp) 2 ' 
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Figure 1: The particle-hole excitation spectra calculated for L = 98, /V^ = Nth = 25,Nt i = 48, and the coupling c = 1.0 
(LEFT) and c = 10 (RIGHT). 



where kf, is the momentum of the hole and k p represent the momentum of quasi-particles. We 
use the same notations Q, D and B to present the new integral limits. 

Now, we calculate the excitation energy E ex = E—Egs , where Ecs is the ground state energy. 
Follow the methods proposed by Takahashi I54l.l55l.l5ql . E/?ler and Korepin 14911 . we define the 
differences of the densities of states between the ground state and the excited state as 

o-m = L[p(k)-p cs (k)lcr 2 (A) = L[pM) -Pic S U)lcr 3 (A m ) = L[p m (A w ) - p™(A w )], (64) 

where pes, Pi.cs and are the corresponding densities at the ground state. The corrections to 
the densities are 



1 



tp\ = -5{k - k h ), ip 2 = - 
Define the bare energies as 



2c 



nc 2 + 4(A - k p ) 2 ' 



<P3 = 0. 



(65) 







f 1 - a f 


o ) 


-i 


' kl 1 


S2 




— a® 


— fij 







«3 , 




, -af 


1 J 




{ o J 



(66) 



where the integral operator a^(x) satisfies 
1 2nc 



a A n (x) 



n n 2 c 2 + Ax 2 ' 
We obtain the excitation energy as 



r 1 2nc 
J_ A n n 2 c 2 + 4(x 



4(x-y) 



-f(y)dy. 



3 pa a 
2 I s a {p)(p a 



dfi = e(k p ) - e(k h ), 



(67) 



(68) 
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Figure 2: The dispersion relations of the elementary excitation: charge (LEFT) and hole (RIGHT). 

where e(k) is the dressed energy 



n J-b c 2 + 4(k - 



(A)dA 



(69) 



1 C Q 2nce(k)dk 1 C B , , , 

- D 2ncef ) (^)^ 1 ' 



1 2ncej 
I f * 

7!" J-B < 



4(^-^ 1 )) 2 ' 



(70) 



The dress energies at zero temperature are calculated from the thermodynamic Limit. The de- 
tailed derivations and the further explanations can be found in the next section. 

5.2. Spin wave excitation 

The second class excitation is flipping one "spin", which means add two holes in the distri- 
bution of the quantum integer series {J„ '}. Then the quantum numbers change from integer to 
half-odd-integer or vice versa. We denote the spectral parameters corresponding to the missing 
integers in the {ji^} by A\ and A h T In the thermodynamic limit we obtain following coupled 
integral equations for a state with two holes 



m=— I f 2c Pl (A)dA 
P{ > 2tt n J_ B c 2 +4(k-A) 2 ' 

1 f 6 2cp(k)dk 1 f 



B cpi(A')dA' 



_ Q c 2 + 4(A - k) 2 7i J_ B c 2 + (A- A') 2 
D 2cp^(A (1 >)dA (1 > 



i r 

n J-D 



+ 4(A-AW) 2 L . , 

7=1 



(1) (1) = 1 f B 2c Pl (A)dA 
P {A > „J B c 2 + 4(A^-A) 2 - U ' 
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Figure 3: The spinon-spinon excitation spectra calculated for L = 98, N/-, = Nf 2 = 25, = 48, and the coupling 
c = 1.0 (LEFT) and c = 10 (RIGHT). 



The differences of the densities of states and the bare energies are defined as before. The 
corrections to the densities at present case are 



0, ( P2 = -Y J S{A-A h j ), V3 =0. 



(73) 



./=! 



We obtain the excitation energy as 



E-E, 



GS 



3 s*a a 

Z \ Ea 



-e(4) - e(4). 



(74) 



This excitation is gapless. The total 5 = 1 and the total S z = 1, thus it is a spin triplet excitation. 

Another spin excitation is the spin singlet excitation. This case means digging two holes in 
the quantum number sequence {Ja 1 } and constructing one /i-string of length 2 (with the quantum 

(2) (2) 

number J a ). Thus in the rapidities {A}, two of them form a 2-string (A. ± ic/2) and the rest are 

(2) 

real. The string center A\ is determined by 

N N h +N fl -2 N bi 

^(Af-k^lnjf^ J] Hl/2 (^-Aa + ^ Sl (^-fl, (75) 

j=\ fi=\ y=l 

If the system-size keeps finite, the Bethe ansatz equations can be solved numerically. The exci- 
tation spectra and dispersion relation are shown in Figs. [3]and|4j respectively. If the system-size 
tends to infinity, we obtain the following coupled integral equations for a state with two holes 
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Figure 4: The dispersion relation of the spinon excitation. 



and one 2-string 

p(k) = — + - — 
2n n J_ B c l + 



2cpi(A)dA 1 



_ B c- -r 4(k - A) 2 nL c 2 + (k- A) 2 ' 

_ i r Q 2c P (k)dk i r B C pi(A')dA' i r u 

~ n J_ e c 2 +4(A-k) 2 ~ n J_ B c 2 +(A- A') 2 + n J_ D 

2c 



7=1 

9 (A ^ nJ_ B c 2 + 4(, 



° 2cp {x \A^)dA^ 
c 2 +4(A-AW) 2 



+ 



6c 



c 2 + 4(A - Af) 2 9c 2 + 4(A - if) 2 , 
1 c 



■40V -A) 2 nL c 2 + {A f _ A) 2- 
The corrections to the densities of states at present case are 



(76) 



1 

<Pi = - 



1 



n c 2 + (k - A) 2 ' 



¥3 



ip 2 



*c 2 + (A\ z> - A) 2 
2c 



6c 



c 2 + 4(A - Af } ) 2 9c 2 + 4(A - A ( . 2) ) 2 , 



(77) 



The excitation energy is E ex = £^=i f_° s a {p)(p a dp. In order to calculate the excitation energy, 
we recall the definition (T70b of the dressed energy e n (A) of /t-string with the length n. It is easy 
to proved that eiiA) - for the present case, then we obtain the excitation energy as 



E ex = -e(A\) - e(4). 



(78) 



We see that the excitation energy (T78T > is the same as that of the excitation with only two holes 
and no strings. 
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Figure 5: The isospinon-isospinon excitation spectra calculated for L = 98, Nf = 50, N\, = 48, and the coupling c = 1.0 
(LEFT) and c = 10 (RIGHT). 



5.3. Isospin excitation 

The third excitation is replacing p holes in the ground state distribution of the quantum num- 
ber series {J^} by a A^-string with the length p, Aj^ = k + ic(p + 1 - 2j)/2, where k is a real 
number and j = 1, • • • , p. The ^-string solution describes an excitation of p particles bound state. 
This excitation is spinless. The quantum numbers / and J a are the same as that at the ground 
state, while the quantum number jumps from half-odd integer to integer. The allowed range 
of integer for the /l (2) -string is |/ (2) | < \(Nb\ - 2p). The center k of the /l (2) -string is determined 
by 



A'/., 



2^ = £ S/)/2 (k-4>>). 



(79) 



7=1 



Please see the next section for further explanations. The densities of states at this excitation are 

rw = _L I f B 2c P^ dA 
P(> 2n nJ_ B c 2 +4(k-A) 2 ' 



1 r e 2cp(k)dk 1 r B c Pl (A')dA' 1 r D 2cp (l \A m )dA m 
~ n J_n c 2 +4(A-k) 2 ~ n J_ B c 2 + (A - A') 2 + n J_ D c 2 + 4(A — A^) 2 ' 



_ e c 2 + MA - k) 2 
B n, — {A)dA 



n J-B c 1 + 4(, 



(A- A') 2 
1 2pc 



_ B c 2 + 4(AW - A) 2 nL p 2 c 2 + 4(A<D - k) 2 ' 
Thus the corrections to the densities are 

1 2pc 

n p 2 c 2 + 4{A m - k) 2 

Using the similar method, we obtain the excitation energy as 



n J-B c 2 + 



B c 2 + 4(A m - k) 2 
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Figure 6: The dispersion relation of the elementary excitation: isospinon. 



Another isospin excitation is the - /l (2) -string excitation. The quantum number series {/,-}, 
{J a } and {/j } are filled symmetrically around the zero in this excitation. In the rapidities 
two of them form the - /l (2) -string, A (l) = /l (2) + ic/2 and the rest are real. The allowed range 
of integer for the A {1) - ^ (2) -string is \J^ } \ < j(N hl - 1). The center of the A m - /l (2) -string k = A {2) 
is determined by 

N tl +Nfy N hl -l 

2njV= Yj Ei(k-Ao)- ^ H,/^*-^). (83) 

a-l y=l 

The Bethe ansatz equations and energy spectrum for the finite system-size case can be solved 
numerically. From that, we obtain the isospinon-isospinon excitation spectra and corresponding 
dispersion relation, which are shown in Figs. [5] and [6] respectively. Comparing Figs. |4]and|6] 
we see that the excitation spectra of spinon is linear while that of isospinon is quadratic for the 
small momentum. 

If the system-size tends to infinity, we obtain the densities of states at this excitation as 



m = L I f 2c Pl (A)dA 
P( > 2k tt J_ B c 2 +4(k-A) 2 ' 



c Pl (A')dA' 



i r 2c P (k)dk 1 r 

P n J_ Q c 2 +4(A-k) 2 n J_ B c 2 + (A - A') 2 

1 r D 2cpV0V)dAW 1 c 
+ n J_ D c 2 + 4(A - A«) 2 + nLc 2 + (A- k) 2 ' 



p (i) U a)) = I f ^ Pl (A)dA _ J_ 2c 

n J- B c 2 +4(A^-A) 2 ttLc 2 + 4(A^-k) 2 ( ' 

Thus the corrections to the densities are 

1 c 1 2c 

<p x = 0, <p 2 = — 2 '. ~7~\ ^> Pi = 2 , At-, ^2" (85) 

n c 1 + (A- k) n c z + 4(A - k) z 
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Using the similar method, we obtain the excitation energy as 

£ «W = - f T—n f 2 777m ^{A m )dA m = 4\k). (86) 

k J-D c 2 +(A- k) 2 n y B c 2 + 40V - k) 2 2 

6. Thermodynamics of the system 

The finite-temperature properties of the system can be studied based on the solutions of the 
Bethe ansatz equations in the thermodynamic limit. The solutions of the Bethe ansatz equations 
are a little bit complicated. Besides real solutions, the Bethe ansatz equations also have complex 
solutions. Generally, the complex solutions are determined by the poles or zeros of the Bethe 
ansatz equations in the thermodynamic limit. Which gives us a hint to determine the structures 
of the solutions of the Bethe ansatz equations. In the following, we only consider the case c > 0, 



because the bosons with attractive interactions (c < 0) do not have the thermodynamics [54]. 
After some analysis, we find that the structures of the solutions of the Bethe ansatz equations 
d49l i - d52l > are: (1) All the momentums kj are real. (2) The rapidities A form the m-string, 
A™ 1 = A[" + (m + 1 - 2j)ic/2 + o(e~ SL ), where X£ are real and j = 1, ■ • • , m. (3) Some rapidities 
A (D _ are real. (4) Some rapidities form the A {1) -/ 2) -strings, 4" = A (V)S + ic/2 + o( e - SL ) 
and Af = where A (1)s are real. (5) Some rapidities /t (2) are real, which are the real part of 
the A m - /l (2, -strings. (6) Some A {2) form the p-strings, Af )pl = Af )p + (p + 1 - 2Z)/c/2 + o( e - SL ), 
where A^ 17 are real and / = 1, • • ■ , p. 

At the temperature T, the system (1) arrive at the thermal equilibrium. We denote the density 
of momentum k by p(k), the density of /l-strings with length n by p„{A), the density of real rapidity 
A m by pf\A m ), the density of A (1) - /l (2) -strings by p 2 V) {A m ) and the density of /l (2) -strings by 
pf(A {2) ). Meanwhile, the notations p h (k), pJ(A), p\ 1)h (A m ), p 2 m (A m ) and pf'\A {2) ) represents 
the densities of corresponding holes. From the Bethe ansatz equations d49l > - d52l . these densities 
should satisfy 



j OO 

p(k) + p\k) = — + Yj a n * Pn(k), 
n=l 

oo 

p„(A) + p h n (A) = a„* p(A) - ^ A„j„ * p m (A) + (A nA - S n> i) * p ( 2 {A) + a„ * pf\A), 

m=l 

oo oo 

p[ l \A^) + p<V) = |>„ *p„(i«) - 2>* *Pf " *P2V°), 

n=\ p=\ 

oo 

pf + p^(A^) = £(A„, - * Pn (A {l) ) ~ar* p^) - a 2 * p?(A<\ 



pf(A^) + pf\A®) = -Y J A M *pf{A^) + a p *p^l (87) 

9=1 

where the integral operators a n (x) and A nm (x) are 
1 2nc 

a„(x) = —z-z —-, floW = o(X), 

n n L c L + 4x z 

A (x)=i a l"- m l w + 2a \n->M+2(x) + ■■■ + 2a n+m - 2 (x) + a n+m {x) if n + m, 
1 2«2(x) + 2ai\(x) + • ■ • + 2a% n -i(x) + ci2 n {x), if n — m. 
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The convolution is denned as a * b(x) = J a(x - y)b{y)dy. 

The Gibbs free energy at a given temperature T, chemical potential A and external magnetic 
field h is given by 

G(T,A,h) = E-AN-h(M f + M b )-TS. (89) 

The densities of energy and particles number are E/L = J Q Q k 2 p(k)dk, N/L — J Q Q p(k)dk. The 
magnetization for fermions and that for bosons are 

OO 

M f = N fl - N f2 = -M\ l) - 2M ( 2 l) + 2 nM„ - N, (90) 

n=\ 

CO 

M h = g(N bl - N h2 ) =2gJ] pMf - £<\ (91) 
where g is the Landau g-factor. The entropy of the system is S/L = j [(p + p h ) ln(p + p h ) 

-plnp - p h lnp h ] dk + Z: =1 / [(Pn + Pi) MPn + A) ~ Pn mp„ - p h „ kip*] dX + / [(p^ + pf) 

ln(pf + p?») - pf > Inpf - p «» Inp^j ^ + J + ^ ln(p (D + ^ _ p U> lnp (D_ 
/£» Inp^ ^(0 + Z- , / [<p«> + p«») ln(p< 2 > + pf ) - p< 2 > lnp< 2 > - p< 2 >* lnp«»] ^. 

For convenience, we introduce following notations £ = p h /p, a n = p h Jp n , Pi = P^ h /p^\ 
/?2 = P2 )h Ip^ an d 7p - P ( p )h /pf^- At the thermodynamic equilibrium state, the free energy must 
be minimized. From the variation of free energy is zero, we obtain following coupled non-linear 
integration equations 

k 2 - A + h A , 
lnf(*) = }]a n *\n(l+a- n l (k)), (92) 

lnar B (A) = --fr + Yj Am '" * ln(1 + ""W) -«»* ln (l+ T'W)) 

m=l 

-a„ * ln(l +P\ l (A)) - (A n>1 - <*„,!) * ln(l + j8^(/l)), (93) 



ln/W (1) ) - - X a n * ln(l + a „-HA (1) )) 



+«! * ln(l + /? 2 V ) )) - J] * ln d + y~ P \^)\ (94) 



ln^ 2 (^ (1) ) = y- Y/An,i - 8 n ,d * ln(l + «; 1 (^ (1) )) 



9=1 

where * denotes the convolution. 



+a 2 * ln(l + /3 2 1 (i (1) )) + ai * ln(l + ft" 1 ^)), (95) 

r 



ln r ^ 2 >) = + £ A 9>p * ln(l + y-\A^)) + a p * ln(l + P\\A^)\ (96) 
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7. Special limits 



In principle, the thermodynamic Bethe ansatz equations can not be solved analytically. One 
has to use the numerical simulations or the approximate methods. However, at some special limit 
cases, the thermodynamic Bethe ansatz equations can be solved exactly. 



J 2 ) 



7.1. Zero temperature limit 

We first consider the zero temperature limit. The dressed energies e, e„, ej !) , and ej/' are 

defined as £ = exp(e/T), a„ = exp(e„/r), /3\ = exp^ 1 '/?"), P2 = exp(e^/7j and y p — 

exp(e p 2) /T^. Substituting them into the thermodynamic Bethe ansatz Eqs. d92b - d%l ) and letting 
the temperature tends to zero, we obtain the analytic formulas for the dressed energies 



e(k) = k 2 + 



n J- 



2c£ n (A)dA 



_ B c 2 + 4(k-A) 2 ' 



e„(X) 



If 



Q 2nce(k)dk 
_ Q c 2 + 4(A - k) 2 



If 

n J-i 



A hn (A - A')ei(A')dA' 



1 r u 2nce\ 
n J-D n 2 c 2 + 

I f B 

n J-b 



1 C D 2nce ( '\A (1) )dA ( 



2ce\ (A)dA 



4V 0) ) = - f 

7!" J-B 

efW = -I f 
n J-L 



_ B c 2 +4{AW-A) 2, 

cei(A)dA 1 C D 2ce\ l) (A)dA 



c 2 + (A™ - T) 2 



-If 

n J_ D c 2 + 40V - A) 2 ' 



2pce^{A^)dA^ 



J- B c 2 + 4{AW - A^) 2 ' 
These equations have been used in calculating the low-lying excitation energies. 



(97) 



(98) 
(99) 
(100) 
(101) 



7.2. High temperature limit 

If the temperature T tends to infinity, the rapidities a n ,/3 \,f$2 and y p become constants and 
the thermodynamic Bethe ansatz equations become a set of coupled algebraic equations. The 
thermodynamic Bethe ansatz Eqs. d93l and ( t96b reads 



In ar = - In 

2 



1 + a 2 



In ct2 — — In 

2 



(1 + ari)(l + ar 3 ) 



1+/L- 



lnar„ = - ln[(l +a„_i)(l + ar„+i)], n > 3, 



ln Tl = iln[(l+ r2 )(l+ y 87 1 )], hy* = ^lnf(l +r P -i)(l +7 P+ i)], P>2. 



(102) 



The solutions of Eq. ( 11021 ) are 



ff„=g z («)-l, n>2, ai =g(2)/f(0), 
y p =f(p)-l, p>l, A = l/[/ 2 (0)-l], 



(103) 
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where g(n) = (a n - a~ n )/(a - a' 1 ), f(p) = (bd p - b~ l d- p )/(d - cT 1 ) and the parameters a,b,d 
are determined by 

k 2 - A + h vV „ i 
Inf = j. 2u ln(1 +a ' ) ' 

n=\ 

ln(l + ai) = -— + 2 £ ln(l + a„>) - ln(l + f 1 ) - ln(l + P^) - ln(l + /32 1 ), 

n=l 

/'I "W 00 00 

InjSi = ^-^-£ln(l +0 + ln(l+/3 2 1 )-^ lnd+y; 1 ), 

n=\ p=\ 

lnjfc = Y~ 2 Yj ln(1 + a " l) + ln(1 + a i 1} + ln(1 1} + ln(1 +/? 2 1} ' 

n=l 

ln(l + 71) = + 2 £ ln(l + y- x ) + ln(l + y 6r 1 ), 

lim ln^ = ^ j^lnr,^ (1Q4) 

7.3. Weak coupling limit 

If the coupling constant c tends to zero, the thermodynamic Bethe ansatz equations can be 
simplified as 

lnf(*) - - iln[(l + or, (*))(! + yS2 1 «)]> 



r 2 

^1 



lnor!^) = 1 [ln(l + tt2 (A)) - ln(l + , 

lna 2 (A) = \ [ln(l + + ln(l + a 3 (A)) - ln(l + y 8 2 1 (^))] , 

lnor„(/t) = iln[(l +a„_i(A))(l +a„ + i(i))], n > 3, 

lnA(^ (1) ) = -^ln[(l +ri(i (1) ))d +r 1 W (1) ))d + 

ln^ 2 (^ (1) ) = -ln[a 1 (/ 1) )(l 

lnyiU (2) ) = ^ [ln(l + y 2 (A & )) - ln(l h-/^ 2 '))] - 

ln 7p (A (2) ) = 1 ln[(l + y p -M i2) m + J P+ M (2) ))l P>2, 
,. ln«„(fc) 2* lny P (^ (2) ) 2g/j 

lim = — , lim = . (105) 

n—><x> n T p—><x> p T 

7.4. Strong coupling limit 

If the coupling constant c tends to infinity, the momentum rapidities are completely decoupled 
with other rapidities. Meanwhile, the solutions for o-„,/?i,/?2 and y p turn into constants and are 
independent of the variables A, A m and A {2) . The thermodynamic Bethe ansatz equations can be 
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solved analytically in this case. The solutions of Eqs. (l93l and d%b still take the form of ( 11031 ). 
where the boundary conditions (11041 ) become ct\ — fi^ and ln/?i = -|ln(l + y\). Then the 
analytic solutions of thermodynamic Bethe ansatz equations are 



n \2 cosh - — h 2 cosh — 
T T 



( sinh f 



cosh J 
cosh y 



Pi 



4cosh 2 ^-T 



sinh 



T > 



1, n>2, 

2 



sinh y 



1, p > 1. (106) 



8. Conclusions 



In summary, we study the exact solutions of ID mixture of spinor bosons and spinor fermions 
with (5-function interactions. The wave function for the bosonic parties is symmetric while for the 
fermionic parties is anti-symmetric. The global wave function of the system is supersymmetric. 
After obtaining the two-body scattering matrix by using the coordinate Bethe ansatz method, 
we prove that the system is integrable. Then we derive the energy spectrum and the Bethe 
ansatz equations with different gradings by using the graded nest quantum inverse scattering 
or algebraic Bethe ansatz method. Based on the solutions of the Bethe ansatz equations, we 
discuss the ground state properties of the system. We find that if the interactions are repulsive, 
the fermions form spin singlet states and the bosons are polarized, thus the global ground state is 
partial polarized. If the interactions are attractive, the bosons condensed at the zero momentum 
point. We discuss the charge-hole excitation, spin wave excitation, isospin wave excitation and 
corresponding excitation energies very detailed. We also obtain the thermodynamic Bethe ansatz 
equations at finite temperature and find their analytic solutions at some special limit cases. 
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